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Abstract 

The  plasma  flow  switch  (PFS)  utilizes  the 
motion  of  a  high  density  plasma  discharge  to 
accumulate  and  then  release  magnetic  energy  at 
multimegampere,  multimegajoule  levels.1  Energy 
Input  to  a  coaxial  plasma  gun  occurs  in  several 
microseconds  as  the  annular  plasma  discharge  is 
accelerated  axially  along  the  center  conductor. 
Energy  is  released  to  a  load  in  a  few  hundred 
nanoseconds  as  this  plasma  is  expelled  off  the  end 
of  the  center  conductor.  Low  density  plasma 
upstream  of  the  axially-moving  discharge  can  limit 
the  flow  of  electromagnetic  power  to  speeds  related 
to  the  local  Alfven  speed.  The  wave  transit  time 
between  the  load  and  the  magnetic  energy  store  can 
compare  with  the  desired  energy  extraction  time. 
The  effects  of  finite  wave  speed  have  beep 
investigated  using  approximate  acoustic  analyses. 
The  present  paper  utilizes  a  transmission  line  code 
with  varying  time  and  impedance  elements  to  model 
wave  transit  effects  during  PFS  operation  including 
commutation  and  load  dynamics. 

I.  Introduction 

The  PFS  utilizes  a  high  density  plasma 
discharge  to  conduct  current  while  electromagnetic 
energy  is  concentrated  in  an  inductive  store  just 
upstream  of  the  PFS  as  illustrated  schematically  in 
Fig.  la.  The  annular  plasma  discharge  in 
cylindrical  geometry  is  initiated  early  in  time  by 
driving  current  through  an  exploding  foil /wire  array 
In  vacuum.  The  plasma  Is  accelerated  downstream  by 
the  axial  JxB  force  and  magnetic  energy  is 
accumulated  ln~  the  inductor  as  the  current  rises. 
The  plasma  mass  and  plasma  gun  electrode  lengths  are 
chosen  such  that  the  plasma  moves  across  an  opening 
gap  with  a  speed  of  several  cm/us  when  the  inductor 
is  fully  energized.  Multimegampere  level  currents 
can  be  conducted  for  several  microseconds  in  this 
fashion,  inductively  storing  energy  at  the 
multimegajoule  level. 

For  the  purpose  of  this  paper  a  dynamic 
imploding  foil  load  is  assumed  to  be  located  in  the 
opening  gap  as  illustrated  in  Fig.  lb.  As  the  high 
density  plasma  moves  across  the  opening  gap  the 
stored  energy  can  be  released  to  the  load,  however, 
low  density  plasma  trailing  the  current  sheet  limits 
the  flow  of  the  electromagnetic  power  to  speeds 
related  to  the  local  Alfren  speed,  v».  The  wave 
transit  time  between  the  load  and  tne  magnetic 
energy  store  can  be  comparable  with  the  desired 
energy  extraction  time.  Typically  energy  is 
released  on  a  timescale  of  a  few  hundred 
nanoseconds.  This  low  density  plasma  results  from 
current  contacts  at  the  surfaces,  UV-ablation  of 
walls,  and  Rayleigh-Taylor  instability  of  the 
upstream  boundary  of  the  accelerating  current 
sheet.  Because  of  its  low  mass,  the  plasma  rapidly 
expands  radially  inward  into  the  gap  filling  the 
entire  volume  to  the  imploding  foil.  Both  the 
plasma  density,  n,  and  the  magnetic  field,  B,  in  the 
opening  gap  are  then  functions  of  time  with  the  wave 
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Fig.  1  Schematic  of  PFS  a)  run  down  and  b) 
commutation 
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speed  varying  as  v.~  B/n  .  The  effects  of  finite 
wave  speed  during  the  opening  phase  of  the  PFS  have 
been  investigated  in  some  detail  using  approximate 
acoustic  analyses.  For  the  purposes  of  this  paper, 
however,  it  will  be  assumed  that  vA  simply  increases 
linearly  in  time. 

The  present  paper  utilizes  a  transmission  line 
code4  with  varying  time  and  impedance  elements  to 
model  PFS  operation  Including  commutation  and  load 
dynamics.  In  a  transmission  line  code  a  circuit 
element  is  modeled  by  a  characteristic  length  or 
transit  time,  t,  expressed  in  seconds  and  a 

characteristic  Impedance,  Z,  expressed  in  ohms.  The 
inductance  of  the  element  is  given  by  L  =  tZ  and  the 
capacitance  of  the  element  is  given  by  C  =  x/Z. 
Forward  moving  voltage  waves,  Vf,  and  backward 
moving  voltage  waves,  V_,  travel  on  the  transmission 
line  elements  with  amplitudes  determined  by  the 
reflection  and  transmission  coefficients  at  the 
various  junctions  of  elements.  The  voltage  at  any 
point  is  given  by  V  =  V+  +  V  and  current  is  given 
by  I  =  (V+-V_)/Z. 

During  the  run-down  of  the  PFS  as  shown  in  Fig. 
la,  the  impedance  of  the  cylindrical  transmission 
line  is  assumed  to  be  fixed  while  the  length 
increases  as  the  plasma  is  accelerated  such 
that  L(t)  =  x(t)Z.  A  discussion  of  how  to  model 

this  type  of  transmission  line  element  within  the 
context  of  the  code  (described  in  Ref.  4)  is  found 
in  Section  II.  During  the  commutation  phase  the 
physical  dimensions  of  the  opening  gap  do  not 
change,  however,  the  wave  transit  speed  through  the 
low  density  plasma  Is  much  less  than  the  speed  of 
light  and  also  varies  with  time.  This  region  then 
is  modeled  by  a  fixed  Inductance,  L,  and  a  variable 
transit  time,  t(t),  such  that  Z(t)”L/x(t).  A 
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II.  Modeling  PFS  Run-Down 


During  the  run-down  of  the  PFS,  it  is  assumed 
that  the  plasma  can  be  approximated  by  a  thin 
current  carrying  annulus  of  mass,  m,  which  is 
uniformly  accelerated  downstream  by  the  JxB  force. 
The  equation  of  motion  for  the  annulus  i s 
2 

m  =  ZfI2/2c,  (1) 

dt 

where  Zf  is  the  fixed  impedance  of  the  cylindrical 
transmission  line  in  which  the  annulus  is  moving, 
I(t)  is  the  current,  £(t)  =  cx(t)  is  the  distance 
the  annulus  has  traveled  and  c  is  the  speed  of  light 
in  vacuum.  The  instantaneous  inductance  of  the  re¬ 
gion  created  by  the  motion  of  the  PFS  is  L(t)=Z,x(t). 
The  velocity  of  the  moving  annulus,  T  v ( t ) 
=  d£/dt=cdx/dt,  is  related  to  the  changing 
inductance  by  dL/dt  =  Z,dx/dt  =  (Z,/c)  v(t).  The 
kinetic  energy  impartedTto  the  plasma  annulus  then 
is  related  to  the  changing  inductance  through 


-  f1  dL  T2Hf 

"  -*2  dt  1  dt’ 


(2) 


r  (timesteps) 


Fig.  2  Illustration  of  redistribution  where  solid 
(dashed)  lines  divide  portion  of  F  in  each  segment 
before  (after)  change  in  length  from  x  to  x. 

transmission  line  junctions  at  the  two  ends  of  the 
element  are  unaltered  because  Z^  is  fixed  in  this 
case. 


The  rate  of  increase  of  electromagnetic  energy 
stored  in  the  transmission  line  element  then  is 
reduced  by  the  rate  of  increase  of  the  kinetic 
energy  of  the  plasma. 

In  order  to  model  the  run-down  of  the  PFS,  a 
transmission  line  element  of  fixed  impedance,  Zf, 
and  varying  length  is  used.  Since  numerically  the 
element  is  an  integral  number  of  timesteps  or 
segments  long  at  any  gven  instant,  the 
length  x  (t)  of  the  element  changes  in  discrete 
steps  inpaccordance  with  Eq.  (1).  Between  timesteps 
the  length  of  the  element  is  instantaneously 
adjusted  such  that  the  charge,  Q  =  C  V,  and 
flux,  <|>  =  LI,  in  the  element  are  conserved.  The 
field  energy  is  not  conserved  in  this  case,  because 
kinetic  energy  is  imparted  to  the  plasma.  When  the 
element  is  N  timesteps  long  (i.e.  x  (t)  =  NAx),  the 

charge  and  flux  in  the  j^*1  segment  of  the  element  at 
the  i*"*1  timestep  are 

v  }  +  v  ! 

Q-=  (—^T - -"J)Ax,  j  =  1, . . .  ,N  (3) 

3  Af 

+]-  (v+j  -  V_J)Ax,  j  =  1 . N  (4) 

where  Ax  is  the  timestep.  If  between  timestep  i  and 
timestep  i  +  1,  the  length  of  the  element  is 

increased  from  NAx  to  MAx,  then  charge  and  flux  are 
distributed  in  the  M  segments  by  maintaining  the 
shape  of  the  distribution  of  Q  and  <|>  along  the 
length  of  the  element  as  illustrated  in  Fig.  2. 
Before  the  voltage  waves  are  advanced  at  timestep  i 
+  1,  the  voltage  waves  on  the  element  are  adjusted 


according  to 

,  i  .  Q]zf 

+j  2  Ax 

A1’ 

<P  .• 

+  2Ax  ’ 

j  =  1 . M 

(5) 

-j  2  Ax 

2Ax  * 

j  =  1 J  •  •  •  jM 

(6) 

where  the  tilde  over  the  variable  indicate  they  are 
the  new  adjusted  variables.  Note  that  reflection 
and  transmission  coefficients  associated  with  the 


As  an  example  consider  the  simple  circuit  shown 
in  Fig.  3.  The  first  element  represents  a  capacitor 
with  C  =  xc/Zc  =  10  uF,  the  second  element  is  an 

inductor  with  Lc  =  Z  x  =  40  nH,  and  the  third 
s  s  s 

element  represents  a  PFS  with  an  initial  inductance 
of  Lp(o)  =  Zfxp(o)  ,=  0.5  nH.  The  circuit  is 

terminated  on  the  left  by  a  large  resistance,  Rc,  to 
simulate  an  open  circuit  and  terminated  on  the  right 
by  a  short  circuit  simulating  the  PFS  plasma.  The 
quarter  cycle  current  risetime  is  1  us  without 
accounting  for  the  PFS  motion  and  is  indicated  by 
the  solid  curve  in  Fig.  4a.  For  a  PFS  mass  of  m  = 
1.6  mg  and  a  travel  of  10.5  cm  or  14  timesteps 
with  Ax  =  0.025  ns,  the  current  is  given  by  the 
dashed  curve  in  Fig.  4a  and  the  PFS  position  is 

given  by  the  dotted  curve.  Figure  4b  shows  the 

energies  stored  in  the  capacitor  and  in  the  inductor 
[i.e.  Ls  +  L  ( t ) ]  and  the  kinetic  energy  of  the 

PFS.  The  peak  current  decreases  when  a  PFS  is 

included  in  the  circuit  because  of  the  increased 

inductance  and  the  energy  lost  to  kinetic  energy  of 
the  PFS  mass. 

III.  Modeling  PFS  Current  Commutation 

As  the  high  density,  high  speed  plasma  moves 

across  the  opening  gap,  current  is  carried  by  a 

trailing  low  density  plasma  as  shown  schematically 
in  Fig.  lb.  The  characteristic  speed  of  the  current 
propagation  through  such  low  density  plasma  is  the 

Alfven  speed  which  is  much  less  than  the  speed  of 

light  in  vacuum.  Although  the  inductance  of  the  gap 
region,  L,,  is  small,  the  transit  time,  x^,can  be 
long,  ana  hence  the  effective  impedance,  Z1?  is 

small.  For  an  imploding  foil  load  this  initial 

transit  time  can  be  comparable  with  the  implosion 

time.  However  as  the  magnetic  field  increases  and 
plasma  density  decreases,  the  transit  time  can 

significantly  decrease,  allowing  the  stored  energy 
to  be  rapidly  released. 

In  order  to  model  this  type  of  transmission 
line  element,  it  is  assumed  that  is  fixed  and 
therefore  Z^(t)  =L,/x1(t).  The  transit  time  can 
be  modeled  by  using  appropriate  acoustic  analysis2. 
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Fig.  4  Plot  of  a)  current  with  no  PFS  motion 
(solid),  with  PFS  motion  (da'shed)  and  PFS  position 
(dotted)  and  b)  energy  distribution  with  PFS  motion. 

however,  for  simplicity  it  is  assumed  here 
that  T.(t)  decreases  linearly  in  time  from  100  ns  to 
5  ns  in  ~  200  ns.  Again  numerically  the  element  is 
an  integral  number  of  timesteps  or  segments  long  at 
any  instant  and  thus  x^(t)  changes  in  discrete 
steps.  The  length  is  adjusted  instantaneously 
between  timesteps.  Unlike  the  case  previously 

discussed  where  kinetic  energy  is  extracted  from  the 
field  energy,  here  the  energy  in  the  element  is 
conserved  as  Zj(t)  and  ^(t)  vary.  When  the 
element  is  N  timesteps  long  (i.e.  x,  (t )  =  NAx),  the 
forward  flowing  energy,  E+,  and  the  Backward  flowing 
energy,  E_,  in  the  jCT  segment  of  the  element  at  the 
itn  timestep  are 
.  2 

E+j  =  (V+j)  Ax/Z^t.),  j  =  1 . N  (7) 

.  2 

E.j  =  (V_j)  Ax/Zi^),  j  =  1 . N  (8) 

where  again  Ax  is  the  timestep.  If  between  timestep 
i  and  1+1,  the  length  of  the  element  is  decreased 
from  NAx  to  MAx  ,  then  the  forward  flowing  energy 
and  backward  flowing  energy  is  distributed  in  the  M 
segments  in  the  same  manner  as  discussed 
previously.  Before  the  voltage  waves  are  advanced 
at  timestep  1+1,  the  voltage  waves  on  the  element 
are  adjusted  according  to 

vj  =  (W^/Ax)1'2,  j  =  1 . N  (9) 

V_j  =  (Z^t.^j/Ax)1/2,  j  =  1 . M  (10) 

where  Z^(t^)  =  NZj(t^)/M  and  again  the  tilde 
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Fig.  5  Circuit  for  PFS  current  commutation. 

indicates  the  adjusted  variables.  In  this  case  care 
must  be  taken  to  preserve  the  sign  of  the  voltage 
waves.  Since  Z-,  changes,  the  reflection  and 
transmission  coefficients  associated  with  the 
transmission  line  junctions  at  the  two  ends  of  the 
element  must  also  be  adjusted  accordingly. 


As  an  example  consider  the  simple  circuit  shown 
in  Fig.  5.  The  first  element  represents  an  inductor 
with  Ln=  Z  x  =  20.  nH,  the  second  element  represents 
the  PFS  wixh  L^Z,  (o)x1(o)  =  2  nH,  and  the  third 
element  represents  the  imploding  load  with 
l2(o)  =  Z„(o)x„=  0.02  nH  and  mass  per  unit 
length  m/x  =  170  mg/cm.  Since  the  intent  here  is 
to  study  transit  time  effects  on  energy  transfer, 
all  three  elements  are  initialized  carrying  10  MA  of 
current.  Most  of  the  energy,  however,  resides  in 
the  storage  inductor  LQ.  The  load  is  modeled  by  a 
simple  one  timestep  long  element  which  has  a  varying 
impedance  to  simulate  the  varying  inductance  of  the 
element55.  The  instantaneous  inductance  (and 
therefore  impedance)  of  the  load  is  calculated  from 
L  =  2x10" '  a  Xn  (r  / r ( t ) )  and  adjusted  between 


timesteps  using  the  simple  slug  model 


m  d2r 
1  dt^ 
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to  determine  the  load  dynamics.  The  voltage  waves 
are  adjusted  in  order  to  conserve  charge  and  flux 
with  energy  extracted  in  the  form  of  kinetic  energy 
of  the  implosion  so  that 


vj=  j  L V^(Z2/Z2+  1)  +  vjCZg/Z 2-  1)]  (12) 


Vl=  j  [vi(Z2/Z2+  1)  +  1)]  (13) 

Because  12  changes,  the  reflection  and  transmission 
coefficients  are  also  adjusted.  In  principle,  the 
model  described  in  the  previous  section  could  be 
used  to  more  accurately  describe  the  imploding 
load.  However,  the  transit  time  of  the  load  can  be 
short  compared  with  that  of  the  PFS  so  the  use  of 
the  numerically  simpler  model  is  warranted  for  the 
purposes  of  this  study. 


Three  cases  have  been  run  for  comparison.  In 
case  (a)  the  transit  time  of  the  PFS,  x,  is  held 
fixed  at  100  ns.  In  this  case  the  energy  stored  in 
the  inductor  LQ  is  prevented  from  being  rapidly 
extracted  to  implode  the  load.  In  case  (b)  the 
transit  time  of  the  PFS  is  held  fixed  at  5  ns.  Here 
the  stored  energy  can  be  rapidly  released  to  drive 
the  load.  In  case  (c)  x^(t)  decreases  linearly  from 
100  ns  to  5  ns  in  ~  200  ns.  This  models  the  actual 
situation  for  a  PFS.  Results  from  these  three  runs 
are  compared  in  Figs.  6,7,8  and  9.  In  Fig.  6 
x  (t)  is  plotted  for  the  three  cases.  In  Fig.  7  the 
load  current  is  plotted,  and  in  Figs.  8  and  9  the 
total  load  energy  and  radius  of  the  imploding  load 
are  plotted.  The  code  runs  terminate,  and  therefore 
the  curves  end,  when  the  imploding  load  reaches  one 
tenth  of  its  original  radius.  Because  of  the 
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proximity  of  the  load  to  the  large  storage  inductor 
in  case  (b),  the  load  current  is  maintained  at  a 
high  value  and  energy  is  readily  extracted  from  the 
store.  In  case  (a),  however,  the  current  drops  as 
energy  is  expended  to  drive  the  load  inward. 
Because  of  the  100  ns  transit  time,  energy  is 
supplied  to  the  load  very  slowly.  In  case  (c), 
which  models  the  PFS,  the  current  initially  drops  as 
in  case  (a)  when  t.~  100  ns.  As  t,  decreases, 
however,  the  current  Degins  to  increase1 and  energy 
is  rapidly  released  to  the  load  as  indicated  in  Fig. 
8.  The  current  actually  overshoots  the  current  in 
case  (b)  because  of  LC  ringing  in  the  simple  circuit 
chosen  to  model  the  problem.  The  finite  capacitance 
of  the  element  modeling  the  PFS  allows  charge  to 
accumulate  as  voltage  is  generated  by  the  imploding 

-1/2 

load.  At  early  time  u>  ~  (LqCj)  is  small  but,  as 

the  transit  time  of  the  PFS  element  decreases, 
w  Increases.  When  reaches  5  ns,  the  period  of 
the  ringing  has  decreased  to  T  =  2*/“  ~  10' 's  close 
to  what  Is  observed.  More  importantly  energy  is 
released  rapidly  from  the  storage  inductor  to  drive 
the  load  as  illustrated  in  Fig.  9.  The  radius  of 
the  Imploding  load  is  also  shown  in  Fig.  9  for 
comparison  of  the  three  cases.  Initially  the  radius 
in  case  (c)  follows  closely  to  that  in  case  (a).  As 
the  transit  time  decreases,  however,  the  implosion 
speed  rapidly  increases  reaching  nearly  the  same 
speed  as  case  (b)  at  the  end  of  the  Implosion. 

IV.  Summary 

Models  for  both  the  PFS  run-down  and  opening 
have  been  presented.  These  models  were  incorporated 
into  a  transmission  line  code  and  were  used  in 
simple  circuits  to  demonstrate  properties  of  PFS 
operation  including  commutation  and  load  dynamics. 
A  more  accurate  treatment  of  the  changing  transit 
time  of  the  PFS  during  commutation  using  approximate 


(c). 

acoustic  analysis  will  be  the  next  improvement.  In 
the  future  these  models  will  be  used  in  more 
sophisticated  circuits  for  comparison  with 
experimental  results  and  for  designing  new 
experiments. 
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